Dense Point Sets Have Sparse Delaunay Triangulations^ 

or ". . . But Not Too Nasty" 

Jeff Erickson 

University of Illinois, Urbana-Champaign 
jeffe@cs.uiuc.edu 
http:/ /www. cs.uiuc.edu/~jeffe 

Submitted to Discrete & Computational Geometry — February 1, 2008 



Abstract 

The spread of a finite set of points is the ratio between the longest and shortest 
pairwise distances. We prove that the Delaunay triangulation of any set of n points 
in with spread A has complexity O(A^). This bound is tight in the worst case 
for all A = 0{^/n). In particular, the Delaunay triangulation of any dense point set 
has linear complexity. We also generalize this upper bound to regular triangulations 
of fe-ply systems of balls, unions of several dense point sets, and uniform samples of 
smooth surfaces. On the other hand, for any n and A = 0(n), we construct a regular 
triangulation of complexity ^}{nA) whose n vertices have spread A. 



* Portions of this work were done while the author was visiting The Ohio State University. This research was 
partiaUy supported by an Sloan Research FeUowship, by NSF CAREER grant CCR-0093348, and by NSF ITR 
grants DMR-0121695 and CCR-0219594. An extended abstract of this paper was presented at the 13th Annual 
ACM-SIAM Symposium on Discrete Algorithms [53]. See http://www.cs.uiuc.edu/~jeffe/pubs/screw.html for the 
most recent version of this paper. 
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1 Introduction 

Delaunay triangulations and Voronoi diagrams are one of the most thoroughly studied objects in 
computational geometry, with applications to nearest-neighbor searching [4, 31, 36, 65], clustering 
[2, 70, 72, 87], finitc-clcmcnt mesh generation [29, 47, 74, 90], deformable surface modeling [28], 
and surface reconstruction [7, 8, 9, 10, 20, 69]. Many algorithms in these application domains begin 
by constructing the Delaunay triangulation or Voronoi diagram of a set of points in M^. Since 
three-dimensional Delaunay triangulations can have complexity Cl{n^) in the worst case, these 
algorithms have worst-case running time 17 (n^). However, this behavior is almost never observed in 
practice except for highly-contrived inputs. For all practical purposes, three-dimensional Delaunay 
triangulations appear to have linear complexity. 

This frustrating discrepancy between theory and practice motivates our investigation of practical 
geometric constraints that imply low-complexity Delaunay triangulations. Previous research on this 
topic has focused on random point sets under various probability distributions [18, 76, 58, 78, 40, 39, 
52, 59, 62]; well-spaced point sets, which are low-discrepancy samples of Lipschitz density functions 
[29, 74, 79, 80, 92, 93]; and surface samples with various density constraints [11, 12, 52, 59, 62]. 
(We will discuss the connections between these models and our results in Section 2.) Our efforts fall 
under the rubric of realistic input models, which have been primarily studied for inputs consisting 
of polygons or polyhedra [16, 100, 101] or sets of balls [66, 101]. 

This paper investigates the complexity of three-dimensional Delaunay triangulations in terms of 
a geometric parameter called the spread, continuing our work in an earlier paper [52] . The spread 
of a set of points is the ratio between the largest and smallest interpoint distances. Of particular 
interest arc dense point sets in M'^, which have spread 0(n^/'^). Valtr and others [6, 49, 96, 97, 98, 99] 
have established several combinatorial results for dense point sets that improve corresponding 
bounds for arbitrary point sets. For example, a dense point set in has at most 0(n^/^) halving 
planes; the best upper bound known for arbitrary point sets is 0(n^/^) [89]. For other combinatorial 
and algorithmic results related to spread, see [5, 23, 33, 55, 63, 64, 71]. 

In Section 3, we prove that the Delaunay triangulation of any set of points in with spread A 
has complexity O(A^). This upper bound is independent of the number of points in the set. In 
particular, the Delaunay triangulation of any dense point set in has only linear complexity. 
This bound is tight in the worst case for all A = 0{^/r^) and improves our earlier upper bound of 
0(A4) [52]. 

Our upper bound can be extended in several ways. To make the notion of spread less sensitive 
to close pairs, we define the order-k spread A^ to be the ratio of the diameter of the set to the 
radius of the smallest ball containing k points. Our proof almost immediately implies that the 
Delaunay triangulation has complexity O(fe^A^) for any k. Our techniques also generalize fairly 
easily to regular triangulations of disjoint balls whose centers have spread A. With somewhat more 
effort, we show that if a set of points can be decomposed into k subsets, each with spread A, 
then its Delaunay triangulation has spread 0{k^A^). Our results also imply upper bounds on the 
complexity of the Delaunay triangulation of uniform or random samples of surfaces. Finally, our 
combinatorial bounds imply that the standard randomized incremental algorithm [65] constructs 
the Delaunay triangulation any set of points in expected time 0(A^ logn). These and other related 
results are developed in Section 4. 

However, our upper bound does not generalize to arbitrary triangulations. In Section 5, for 
any n and A < n, we construct a regular triangulation, whose n vertices have spread A, whose 
overall complexity is ri(nA). (The defining balls for this triangulation overlap heavily.) This 
worst-case lower bound was already known for Delaunay triangulations for all s/n < A < n [52] . 
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In particular, there is a dense point set in M^, arbitrarily close to a cubical lattice, with a regular 
triangulation of complexity il(n^/^). 

Throughout the paper, we analyze the complexity of three-dimensional Delaunay triangulations 
by counting their edges. Since the link of every vertex in a three-dimensional triangulation is a 
planar graph, Euler's formula implies that any triangulation with n vertices and e edges has at 
most 2e — 2n triangles and e — n tetrahedra. Two points are joined by an edge in the Delaunay 
triangulation of a set S if and only if they lie on a sphere with no points of S in its interior. 

2 Previous and related results 
2.1 Points in Space 

Our results for dense sets compare favorably with three other types of 'realistic' point data: points 
with small integer coordinates, random points, and well-spaced points. Figure 1 illustrates these 
four models. Although the results for these models are quite similar, we emphasize that with one 
exception — integer points with small coordinates are dense — results in each model are formally 
incomparable with results in any other model. 

Unlike our new results, which apply only to points in 3-space, all of these related results have 
been generalized to higher dimensions. We conjecture that the Delaunay triangulation of any 
cZ-dimensional point set with spread A has complexity 0{A'^), but new proof techniques will be 
required to prove this bound. 

Integer points. First, we easily observe that any triangulation of n points in with integer 
coordinates between 1 and A has complexity O(A^), since each tetrahedron has volume at least 1/6. 

(It is an open question whether this bound is tight for all n and A.) In particular, if A = 0(n^/^), 
so that the set is dense, the complexity of the Delaunay triangulation is 0{n). Dense sets obviously 
need not lie on a coarse integer grid; nevertheless, this observation provides some useful intuition 
for our results. 

Random points. Statistical properties of Voronoi diagrams of random points have been studied 
for decades, much longer than any systematic algorithmic development. In the early 1950s, Mei- 
jering [76] proved that for a homogeneous Poisson process in R^, the expected number of Delaunay 
neighbors of any point is 487r^/35-|-2 Ri 15.54; see Miles [78]. This result immediately implies that 
the Delaunay triangulation of a sufficiently dense random periodic point set has linear expected 
complexity. See M0ller [82] for generalizations to arbitrary dimensions and Okabe et al. [86, Chap- 
ter 5] for an extensive survey of statistical properties of random Voronoi diagrams and Delaunay 
triangulations. 

Bentley et al. [14] proved that for n Tiniformly distributed points in the d-dimensional hyper- 
cube, the expected number of points with more than a constant number of Delaunay neighbors 
is 0(n^~^/'^ log n); all such points lie near the boundary of the cube. Extending their technique, 
Bernal [18] proved that the expected complexity of the Delaunay triangulation of n random points 
in the three-dimensional cube is 0{n). Dwyer [40, 39] showed that if a set of n points is generated 
uniformly at random in the unit ball in W^, the Delaunay triangulation has expected complexity 
(^0(c/)j^. particular, each point has dP^'^^ Delaunay neighbors on average. 

Random point sets are not dense, even in expectation. Let 5 be a set of n points, generated 
independently and uniformly from the unit hypercube in W^. A straightforward 'balls and bins' 
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(c) (d) 
Figure 1. Four models of 'realistic' point sets, (a) Small integer, (b) Random, (c) Well-spaced, (d) Dense. 

argument [50, 83] implies that the expected spread of S is 0(n^/^); moreover, with high probabihty,"*^ 
the spread of S is between r2(ra^/'^/ log^ n) and 0(n^/'^"'"^) for any e. 



Well-spaced points. Miller, Talmor, Teng and others [79, 80, 92, 93] have derived several results 
for well-spaced point sets in the context of high-quality mesh generation. A point set S in is 
well-spaced with respect to a 1-Lipschitz spacing function A : M'^ ^ M"'" if, for some fixed constants 
< 6 < 1/2 and < e < 1, the distance from any point x G M"^ to its second nearest neighbor^ 
p G 5 is between 6eX{p) and e\{p). In her thesis, Talmor [92] proves that Delaunay triangulations 
of well-spaced point sets have complexity 0(n); in particular, any point in a well-spaced point set 
has 0(1) Delaunay neighbors. 

Any point set that is well-spaced with respect to a constant spacing function is dense. In general, 
however, the spread of a well-spaced set can be exponentially large; consider the one-dimensional 
well-spaced set {2~* \ 1 < i < n}. On the other hand, dense point sets can contain large gaps 
and thus are not necessarily well-spaced with respect to any Lipschitz spacing function; compare 
Figures 1(c) and 1(d). 

^that is, with probability 1 — n~'^ for some constant c > 

^If a; is a point in S, then x is its own nearest neighbor in S. This is not the definition actually proposed by Miller 
et al., but it is easy to prove that our definition is equivalent to theirs. 
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Talmor's linear upper bound [92] depends exponentially on the spacing parameters S and £, but 
this is largely an artifact of the generality of her results.^ For the three-dimensional case, we can 
derive tighter bounds as follows. Let 5 be a point set in M.^ that is well-spaced with respect to 
some 1-Lipschitz spacing function A. Let p be any point in S, and let q and r be any two Delaunay 
neighbors of p. We immediately have the inequalities \pq\ < mm{eX{p),eX{q)} and \qr\ > 5eX{q). 
Because the spacing function A is Lipschitz, we have X{q) > X(p)— \pq\ > X{p) — eX{q), which implies 
that X{q) > X{p)/{1 + s). Together, these inequalities imply that the set of Delaunay neighbors 
of p has spread at most 2(1 + e)/S = 0(1/5). A packing argument in our earlier paper [52] now 
implies that p has 0(1/6'^) Delaunay neighbors. It follows that the Delaunay triangulation of S 
has complexity 0{n/5'^). (Surprisingly, this bound does not depend on e at all!) 

Finally, in contrast to both random and well-spaced point sets, a single point in a dense set 
in can have 6(A^) = 0(n^/^) Delaunay neighbors in the worst case [52]. Thus, our upper bound 
proof must consider global properties of the Delaunay triangulation. 

2.2 Points on Surfaces 

The complexity of Delaunay triangulations of points on two-dimensional surfaces in space has 
also been studied, largely due to the recent proliferation of Delaunay-based surface reconstruction 
algorithms [7, 8, 9, 10, 20, 69]. Upper and lower bounds range from linear to quadratic, depending 
on exactly how the problem is formulated. Specifically, the results depend on whether the surface 
is considered fixed or variable, whether the surface is smooth or polyhedral, and on the precise 
sampling conditions to be analyzed. We first review some standard terminology. 

Let S be a (7^ surface embedded in M^. A medial hall of S is a ball whose interior is disjoint 
from E and whose boundary touches S at more than one point. The center of a medial ball is 
called a medial point, and the closure of the set of medial points is the medial axis of S. The local 
feature size of a point x G S, denoted lfs(a;), is the the distance from x to the medial axis. Finally, 
a set of points P C S is an e-sample of S if the distance from any surface point x G E to the 
nearest sample point in P is at most £-lfs(a;). This condition imposes a lower bound on the number 
of sample points in any region of the surface. Given an e-sample of an unknown surface S, for 
sufficiently small e, the algorithms cited above provably reconstruct a surface geometrically close 
and topologically equivalent to E. As a first step, each algorithm constructs the Voronoi diagram of 
the sample points; the complexity of this Voronoi diagram is clearly a lower bound on the running 
time of the algorithm. 

Unfortunately, £-samples can have arbitrarily complex Delaunay triangulations due to oversam- 
pling. Specifically, for any surface other than the sphere and any sampling density e > 0, there is 
an e-sample whose Delaunay triangulation has complexity 0(n^), where n is the number of sample 
points [52]. Thus, in order to obtain nontrivial upper bounds, we must also impose an upper bound 
on the density of samples. 

Dey, Funke, and Ramos [37, 54] define"^ a set of points P C S to be a locally uniform sample 
of S if P is well-spaced (in the sense of Miller et al.) with respect to some 1-Lipschitz function 
A : E ^ M+ such that X{x) < Ifs(x) for all a; G E. If P is well-spaced with respect to the local feature 
size function, which is always 1-Lipschitz, we call P a uniform sample of S [52]. Dey et al. [37] 
described an algorithm to reconstruct a surface from a locally uniform sample in 0(n log n) time; 

^Specifically, Talmor first proves that the Delaunay triangulation of any well-spaced point sets in any fixed dimen- 
sions is well-shaped: for every simplex, the ratio of circumradius to shortest edge length is bounded by a constant. 
She then proves that any vertex in any well-shaped triangulation is incident to only a constant number of simplices. 
Again, this is not the definition proposed by Dey et al, but it is easy to show that the definitions are equivalent. 
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Funke and Ramos later showed to how extract a locally uniform e-sample from an arbitrary e- 
sample in 0(n log n) time. Neither of these algorithms constructs the Delaunay triangulation of 
the points. 

In our earlier paper [52], we derived lower bounds on the complexity of Delaunay triangulations 
of uniform samples in terms of the sample measure 



Any uniform e-sample of S contains 0(e^^(S)) points. There are smooth connected surfaces with 
sample measure /x for which the Delaunay triangulation of any uniform e-sample has complexity 



More positive results can be obtained by considering the surface to be fixed and considering 
the asymptotic complexity of the Delaunay triangulation as the number of sample points tends to 
infinity. In this context, hidden constants in the upper bounds depend on geometric parameters 
of the fixed surface, such as the number of facets or their maximum aspect ratio for polyhedral 
surfaces, or the minimum curvature radius or sample measure for curved surfaces. Since the surface 
is fixed, all such parameters are considered constants. 

Golin and Na [59, 60, 61] proved that if n points are chosen uniformly at random on the 
surface of any fixed three-dimensional convex polytope, the expected complexity of their Delaunay 
triangulation is 0{n). Using similar techniques, they recently showed that a random sample of a 
fixed nonconvex polyhedron has Delaunay complexity O(nlog^n) with high probability [62]. In 
fact, their analysis applies to any fixed set of disjoint triangles in R^. 

Attali and Boissonnat [12] recently proved that the Delaunay triangulation of any (e, k) -sample 
of a fixed polyhedral surface has complexity 0{K?n), improving their previous upper bound of 
0{rJ/'^) (for constant k) [11]. A set of points is called an (e, K)-sample of a surface S if every ball 
of radius e whose center lies on S contains at least one and at most k points in P? A simple 
application of Chernoff bounds (see Theorem 4.11) implies that a random sample of n points on a 
fixed surface is an 0(e, 0(log n))-sample with high probability, where e = 0{-\/ (log n)/n). Thus, 
Attali and Boissonnat 's result improves Golin and Na's high-probability bound for random points 
to 0(n log^ n). 

Our new upper bound has a similar corollary. Informally, a uniform sample of any fixed (not 
necessarily polyhedral, smooth, or convex) surface has spread 0{-s/n), so its Delaunay triangulation 

has complexity 0(n^/^). This bound is tight in the worst case; a right circular cylinder with constant 
height and radius has a uniform (e, l)-samplc with Delaunay complexity Jl(n^/^). Similar arguments 
establish upper bounds of 0(«;^n^/^) for (e, /t)-samples and 0(n'^/^ log'^^^ n), with high probability, 
for random samples. We describe these results more formally in Section 4. 

3 Sparse Delaunay Triangulations 

In this section, we prove the main result of the paper. 

Theorem 3.1. The Delaunay triangulation of any finite set of points in with spread A iias 
complexity O ( ) . 

®This definition ignores the local feature size, which is necessary for polyhedral surfaces, since the local feature 
size is zero at any sharp corner. Moreover, for any fixed smooth surface, the minimum local feature size is a constant, 
so any (g, /t)-sample is an 0(e)-sample according to our earlier definition. 




0(/xVlogV)- 
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Our proof is structured as follows. We will implicitly assume that no two points are closer than 
unit distance apart, so that spread is synonymous with diameter. Two sets P and Q are well- 
separated if each set lies inside a ball of radius r, and these two balls are separated by distance 2r. 
Without loss of generality, we assume that the balls containing P and Q are centered at points 
(0, 0, 2r) and (0, 0, — 2r), respectively. Our argument ultimately reduces to counting the number of 
crossing edges — edges in the Delaunay triangulation of P U Q with one endpoint in each set. See 
Figure 2. 




Figure 2. A well-separated pair of sets P(J Q and a crossing edge intersecting a pixel. 
Our proof has four major steps, each presented in its own subsection. 

• We place a grid of O(r^) circular pixels of constant radius e on the plane z = 0, so that every 
crossing edge passes through a pixel. In Section 3.1, we prove that all the crossing edges 
stabbing any single pixel lie within a slab of constant width between two parallel planes. Our 
proof relies on the fact that the edges of a Delaunay triangulation have a consistent depth 
order from any viewpoint [42]. 

• We say that a crossing edge is relaxed if its endpoints lie on an empty sphere of radius 0(r). 
In Section 3.2, we show that at most 0{r) relaxed edges pass through any pixel, using a 
generalization of the 'Swiss cheese' packing argument used to prove our earlier O(A^) upper 
bound [52]. This implies that there are O(r^) relaxed crossing edges overall. 

• In Section 3.3, we show that there are a constant number of conformal (i. e., sphere-preserving) 
transformations that change the spread of P U Q by at most a constant factor, such that 
every crossing edge of P U Q is a relaxed Delaunay edge in at least one image. The proof 
uses a packing argument in a particular subspace of the space of three-dimensional Mobius 
transformations. It follows that PUQ has at most O(r^) crossing edges. 

• Finally, in Section 3.4, we count the Delaunay edges for an arbitrary point set S using an 
octtree-based well-separated pair decomposition [22]. Every edge in the Delaunay triangula- 
tion of S" is a crossing edge of some subset pair in the decomposition. However, not every 
crossing edge is a Delaunay edge; a subset pair contributes a Delaunay edge only if it is close 
to a large empty witness ball. We charge the pair's O(r^) crossing edges to the il(r^) volume 
of this ball. We choose the witness balls so that any unit of volume is charged at most a 
constant number of times, implying the final O(A^) bound. 

3.1 Nearly Concurrent Crossing Edges Are Nearly Coplanar 

The first step in our proof is to show that the crossing edges intersecting any pixel are nearly 
coplanar. To do this, we use an important fact about depth orders of Delaunay triangulations. 
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related to shellings of convex polytopes. 

Let x be a point in M^, called the viewpoint, and let 5 be a set of line segments (or other 
convex objects). A segment s G S is behind another segment t G S with respect to a: if t intersects 
convja;, s}. If the transitive closure of this relation is a partial order, any linear extension is called 
a consistent depth order of S with respect to x. Otherwise, S contains a depth cycle — a sequence 
of segments si, S2, • • • , Sfc such that every segment Si is directly behind its successor Sj+i and s^ is 
directly behind si. De Berg et al. [17] describe an algorithm that either computes a depth order or 
finds a depth cycle for a given set of n segments, in 0{n'^^^^^) time. See de Berg [15] and Chazelle 
et al. [27] for related results. 

We say that three line segments form a screw if they form a depth cycle from some viewpoint. 
See Figure 3(a). 

Lemma 3.2. The edges of any Delaunay triangulation have a consistent depth order from any 
viewpoint. In particular, no three Delaunay edges form a screw. 

Proof: Let x be a point, and let 5" be a sphere with radius r and center c. The power distance 
from X to 5 is |rccp — r^; if x is outside S, this is the square of the distance from x to 5 along 
a line tangent to S. Edelsbrunner [42, 44] proved that a consistent depth order for the simplices 
in any Delaunay triangulation, with respect to any viewpoint x, can be obtained by sorting the 
power distances from x to the (empty) circumspheres of the simplices. (See Section 4.3.) This is 
precisely the order in which the Delaunay tetrahedra are computed by Seidel's shelling convex hull 
algorithm [88]. We can easily extract a consistent depth order for the Delaunay edges from this 
simplex order. □ 

The next lemma describes sufficient (but not necessary) conditions for three pairwise-skew 
segments to form a screw. 




(c) 

Figure 3. (a) A screw, (b) Front view of C, showing viewpoints where Sx appears behind Sz. (c) Every vertex of C sees 
a different depth order, two of which are inconsistent. See Lemma 3.3. 
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Lemma 3.3. Let c be a parallelepiped centered at the origin, and let C = w ■ c for some w > 
2 + \/5 ~ 4.2361. Three line segments, each parallel to a different edge of c, form a screw if they 
all intersect c but none of their endpoints lie inside C. 

Proof: Since any affine image of a screw is also a screw, it suffices to consider the case where c 
and C are concentric axis-ahgned cubes of widths 1 and w, respectively. Let Sx,Sy, Sz be the three 
segments, parallel to the x-, y-, and z-axes, respectively. Any ordered pair of these segments, say 
{sz, Sx), define an unbounded polyhedral region Vxz of viewpoints from which Sx appears behind s^. 
The segment is the only bounded edge of Vxz, and both of its endpoints are outside C. Thus, 
we can determine which vertices of C lie inside Vxz by considering the projection to the xy-plane. 
From Figure 3(b), we observe that if tt; > 2 + \/5, then Vxz contains exactly half of the vertices 
of C, all with the same y-coordinate. A symmetric argument implies that the other four vertices 
lie in Vzx- Similarly, Vxy and Vyx partition the vertices of C by their z-coordinates, and Vyz and 
Vzy partition the vertices of C along their x-coordinates. Thus, each of the eight vertices of C 
sees one of the eight possible depth orders of the three segments. Since only six of these orders 
are consistent, two vertices of C see a depth cycle, implying that the segments form a screw. See 
Figure 3(c). □ 

Recall that a pixel is a circle of radius e in the xy-plane. 

Lemma 3.4. The crossing edges passing through any pixel lie inside a slab of width (20 + 9a/5)£ ~ 
40.1246e between two parallel planes. 

Proof: We will in fact prove a stronger statement. Any two planes hi and /12 whose line of 
intersection lies in the xy-plane define an anchored double wedge, consisting of all the points above hi 
and below /i2 or vice versa. We define the thickness of an anchored double wedge to be the width 
of the two-dimensional slab obtained by intersecting the double wedge with the plane z = r. We 
claim that the crossing edges passing through any pixel lie inside an e-ncighborhood of an anchored 
double wedge with thickness (6 + 3\/5)e- The lemma follows immediately from this claim. 

Without loss of generality, suppose the pixel tt is centered at the origin, and let E denote the 
set of crossing edges passing through tt. Translate each edge in E parallel to the xy-plane so that 
it passes through the origin, and call the resulting set of segments E. We need to show that the 
segments E lie in an anchored double wedge of thickness (6 + 3-\/5)£- Since all these segments pass 
through the origin, it suffices to show that the intersection points between E and the plane z = r 
lie in a two-dimensional strip of width (6 + 3\/5)£ between two parallel lines. The width of a set 
of planar points is determined by only three points, so it suffices to check every triple of segments 
in E. 

Let 61,62,63 be three arbitrary crossing edges in E, and let 61,62,63 be the corresponding 
segments in E. For each i, let pi and qi be the intersection points of with the planes z = r and 

z = —r, respectively, and let Sj C gj be the segment between pi and qi. Finally, let uj be the width 
of the thinnest two-dimensional strip containing the triangle Apip2P3- Observe that A contains a 
circle of radius u;/3. 

Let C denote the scaled Minkowski sum (si + S2 + S3)/3; this is a parallelepiped centered at 
the origin, with edges parallel to the original crossing edges e^. Since each segment Sj fits exactly 
within the slab — r < z < r, so does the cuboid C. The intersection of C with the xy-plane is the 
scaled Minkowski sum (Apip2P3 + Aqiq2q3)/2; this hexagon also contains a circle of radius uo/S. 

Let c be the smallest parallelepiped similar to and concentric with C that contains the pixel tt. 
Since tt is a circle of radius e in the xy-plane, C is at least a factor of a;/3£ larger than c. Lemma 3.2 
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implies that the Delaunay edges ei, 62, 63 cannot form a screw. Thus, by Lemma 3.3, we must have 
uj/3e < 2 + \/5i or equivalently, u; < (6 + 3-v/5)e, as claimed. □ 

3.2 Slabs Contain Few Relaxed Edges 

At this point, we would like to argue that any slab of constant width contains only 0{r) crossing 
edges. Unfortunately, this is not true — a variant of our helix construction [52] implies that a slab 
can contain up to Q{r^) edges, i7(r^) of which can pass through a single, arbitrarily small pixel. 
However, most of these Delaunay edges have extremely large empty circumspheres. 

We say that a crossing edge is relaxed if its endpoints lie on the boundary of an empty ball 
with radius less than 4r, and tense otherwise. In this section, we show that few relaxed edges pass 
through any pixel. Once again, recall that that e denotes the pixel radius. 

Lemma 3.5. If e < 1/16, then for any pixel tt, each point in P is an endpoint of at most one 
relaxed edge passing through tt. 

Proof: Suppose some point p G P is an endpoint of two crossing edges pq and pq', both passing 
through TT, where \pq\ > \pq'\. We immediately have ^qpq' < 2tan~^(e/r) and \qq'\ > 1. See Figure 
4. Thus, the circle through p, q, and q' has radius at least l/(4tan~^(e/r)) « r/4e > 4r. Any 
empty circumsphere of pq must have at least this radius, so it must be tense. □ 




Figure 4. Proof of Lemma 3.5. The radius of the circle is at least 1/(4 tan ^(e/r)), so edge pq' is tense. 

Lemma 3.6. The relaxed edges inside any slab of constant width are incident to at most 0{r) 
points in P. 

Proof: Let o" be a slab of width lo = 0(1) between two parallel planes, and let E be the set of 
points in P incident to any relaxed edge contained in a. To prove that \E\ = 0{r), we use a 
variant of our earlier 'Swiss cheese' packing argument [52]. Intuitively, we take the intersection 
of the bounding sphere of P and the slab a, remove the Delaunay circumspheres of relaxed edges 
in (T, argue that the resulting 'Swiss cheese slice' has small surface area, and then charge a constant 
amount of surface area to each endpoint in E. To formalize this argument, we need to slightly 
expand a and slightly contract the Delaunay balls. 

Let a' be a parallel slab with the same central plane as a, with slightly larger width a; + 1. Let 
denote the sphere of radius r containing P. Let D be the intersection of a' with the sphere of 
radius r + 1 concentric with QP. The volume of D is at most tt{lo + l)(r + 1)^ = O(r^). 

For each point p € E, we define two balls: Bp is the smallest Delaunay ball of some relaxed 
edge pq, and bp is the open ball concentric with Bp but with radius smaller by 1/3. The radius 
of Bp is at most 4r, and the radius of 5^ is at most 4r — 1/3. 
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Figure 5. The Swiss cheese slice E determined by two relaxed crossing edges intersecting a common pixel. The darker 
(blue) portion of the surface is H. The contracted Delaunay balls bp are not shown. 

Finally, we define the 'Swiss cheese shce' S = D \ |JpG_E ^p- each point p £ E, let hp = 
dTi n dhp be the concave surface of the corresponding 'hole' not eaten by any other ball, and let 
H = UpGE ^p- Equivalently, H = dT,\ dD. See Figure 5. 

We claim that the surface area of H is only 0{r). The proof of this claim is elementary but 
tedious; we give the proof separately below. By an argument identical to Lemma 2.6 of [52], the 
ball of unit diameter centered at any point p € E contains at least ^2(1) of this surface area — for 
completeness, we also include this argument below. Since these unit-diameter balls are disjoint, we 
conclude that E contains at most 0{r) points. □ 

Claim 3.6.1. area(F) = 0(r). 

Proof: Let Cp be the common center of Bp and bp, and let ap be the axis line through Cp and normal 
to the planes bounding a. For any point x E hp, let x be its nearest neighbor on the axis Op, and 
let Hp be the union of segments xx over all x G hp. See Figure 6 for a two-dimensional example. 
Finally, let rp = m.inx£hp\xx\. 

The triangle inequality implies that xx and yy have disjoint interiors whenever x ^ y, so 




(1) 



We can bound the volume of each hole H„ as follows: 




(2) 



The intersection bp D da' consists of two parallel disks, the smaller of which has radius Vp. Since 
both cr' and the boundary of Bp contain the endpoints of some crossing edge pq of length at least 2r, 
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the larger of these two disks has radius at least r — lo — 4/3. Again referring to Figure 6, we choose 
two points y, z £ hpCi da' C 6p n da' on the boundary of the larger and smaller disks, respectively, 
so that r„ = \zz\. We can bound this radius as follows: 



\zz\ 



> 



ZCn 



ZCn 



{\zy\ + \ycp 
zy\ - 



ycp\ 



zy\ +2\zy\ 



\yy\' 



|2 I -|2 , I |2 

ycp\ - \yy\ + \ycp\ 



\yy\ 



yy\ -'2\zy 

-|2 

yy\ 



ycp\ 



\yy\ 



2 \zy\ \ycp\ 



\zy\ 



Now substituting the known equations and inequalities 

\yy\ >r — oj — 4/3, \zy\ = u; + 1, 
we obtain the lower bound 



\ycp\ <Ar- 1/3, 



> (r - 1/3)2 - 2{ijj + l)(6r - 1/3) - {uj + if = Q.{r'^ 



(3) 



Finally, combining inequalities (1), (2), and (3) yields an upper bound for the surface area of H: 



1 2r 

area(ii') = ^ area(/ip) < ^ — vol(Fp) = 0(l/r) • ^vol(Fp) = 0{r). 

peE peE peE 



□ 



Claim 3.6.2 (Lemma 2.6 of [52]). Let U be any unit-diameter ball contained in S whose center 
is distance 1/3 from H. Then U contains f^(l) surface area of H. 

Proof: Without loss of generality, assume that U is centered at the origin and that (0,0,1/3) is 
the closest point of H to the origin. Let U' be the open ball of radius 1/3 centered at the origin, 
let V be the open unit ball centered at (0, 0, 5/6), and let W be the cone whose apex is the origin 
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and whose base is the circle dU D dV. See Figure 7. U' lies entirely inside S, and since r > 1, we 
easily observe that V lies entirely outside S. Thus, the surface area of H nW H nU is at least 
the area of the spherical cap dU' Pi W, which is 7r/27. □ 

Together, Lemmas 3.4, 3.5, and 3.6 imply that 0{r) relaxed edges intersect any pixel. Since 
there are O^r"^) pixels, we conclude that there are 0{r^) relaxed edges overall. 

3.3 Tense Edges Are Easy to Relax 

In order to count the tense crossing edges of P U Q, we will show that there are a constant number 

of transformations of space, such that every tense edge is mapped to a relaxed edge at least once. 

' — -1 

A Mobius transformation is a continuous bijection from the extended Euclidean space M = 
M.'^ U {oo} ~ S"^ to itself, such that the image of any sphere is a sphere. (A hyperplane in is a 
sphere through oo in M .) The space of Mobius transformations is generated by inversions. Exam- 
ples include reflections (inversions by hyperplanes), translations (the composition of two parallel 
reflections), dilations (the composition of two concentric inversions), and the well-known stere- 
ographic lifting map from to E>'^ C M'^"'"^ relating d-dimensional Delaunay triangulations to 
{d + l)-dimensional convex hulls [21]: 

X{xi,X2, ...,Xd)= ^^dA) ^ y^^^-^ = (0, 0, . . . , 0, 0). 

xi + x^^ h 1 

Mobius transformations are also the maps induced on the boundary of hyperbolic space S'^^^ by 
hyperbolic isometrics. Two-dimensional Mobius transformations are also called linear fractional 
transformations^ since they can be written as maps on the extended complex plane C = CU {oo} ~ 
§^ of the form z ^ {az + b)/ {cz + d) for some complex numbers a, 6, c, d. 

Mobius transformations are conformal, meaning they locally preserve angles. There are infinitely 
many other two-dimensional conformal maps [85] — in fact, conformal maps are widely used in algo- 
rithms for meshing planar domains [38] and parameterizing surfaces [41, 73] — but Mobius transfor- 
mations are the only conformal maps in dimensions three and higher. Higher-dimensional Mobius 
transformations are described in detail by Beardon [13]; see also Hilbert and Cohn-Vossen [68], 
Thurston [95, 94], or Miller et al. [81]. 

Let S be a sphere in with finite radius (not passing through the point oo), and let vr : ^ 
be a conformal transformation. If '7r(S) is also a finite-radius sphere and the point 7r(oo) lies in the 
interior of n(S), we say that vr everts S. 

Let S he a set of points in C M^, and let p,q,r,s €^ S be the vertices of a Delaunay simplex 
with empty circumsphere S. For any conformal transformation k, the sphere k(E) passes through 
the points k{p), n{q), ^(r), and k{s). This sphere either excludes every other point in n{S), contains 
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every other point in k{S), or is a plane with every other point of k(S) on one side. In other words, 
con'v{K{p), K{q), K{r), k{s)} is either a Delaunay simplex, an anti-Delaunay^ simplex, or a convex 
hull facet of k{S). Thus, ignoring degenerate cases, the abstract simplicial complex consisting of 
Delaunay and anti-Delaunay simplices of any point set, which we call its Delaunay polytope, is 
invariant under conformal transformations. 

In this section, we exploit this conformal invariance to count tense crossing edges. The main 
idea is to find a small collection of conformal maps, such that for any tense edge, at least one of 
the maps transforms it into a relaxed edge, by shrinking (but not everting) its circumsphere. In 
order to apply our earlier arguments to count the transformed edges, we consider only conformal 
maps that map PU Q to another well-separated pair of sets with nearly the same spread. 

Recall that P and Q lie inside balls of radius r centered at (0, 0, 2r) and (0, 0, — 2r), respectively. 
Call these balls QP and Q)Q. We call an orientation-preserving conformal map k a rotary map 
if it preserves these spheres, that is, if k{Q)P) = Q)P and k{Q)Q) = Q)Q. Rotary maps actually 
preserve a continuous one-parameter family of spheres centered on the z-axis, including the points 
p* = (0, 0, \/?>r) and q* = (0, 0, — \/3r) and the plane z = 0. (In the space of spheres [35, 43], this 
family is just the line through Q)P and Q)Q.) See Figure 8 for a two-dimensional example. 




Figure 8. Every rotary map preserves every solid circle and maps every dotted circle to another dotted circle. The bold 
circles are QP and OQ- 

The image of PUQ under any rotary map is clearly well-separated. In order to apply our earlier 
arguments, we also require that these maps do not significantly change the spread. 

Lemma 3.7. For any rotary map k, the closest pair of points in k{P U Q) has distance between 
1/3 and 3. 

Proof: Consider the stereographic lifting map A : M ^ S that takes p* and q* to opposite 
poles of and the plane z = to the equatorial sphere. Any rotary map can be written as 

o po A, where p is a simple rotation about the axis X{p*)X{q*). (Thus, the space of rotary maps 
is isomorphic to SO (3), the group of rigid motions of S^.) 

To make the stereographic lifting map A concrete, we embed and S'^ into M^, as the hyperplane 
X4 = \/3r and the sphere of radius \/3r/2 centered at {0,0,0, VSr /2), respectively. Now A is an 

®The anti-Delaunay triangulation is the dual of the furthest point Voronoi diagram. 
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inversion through the sphere of radius \/3r centered at the origin o = (0,0,0,0): 

3r^(xi,X2,X3,a;4) 
A(a;i,a;2,a;3,a;4) = ^ ' ^ ^ ^ . A(oo) = (0, 0, 0, 0) 

Simple calculations (see Beardon [13, pp. 26-27]) imply that for any points p,q & M , we have 

3r^ \pq\ 



\\{p)X{q)\ 



\po\ \qo\ ■ 



The distance from the origin o to any point in P U Q (in the hyperplane w = \/3r) is between 2r 
and 2\/3r. Thus, for any points p,qEPUQ,we have 



4 - ' - 4 

Simple rotations do not change distances at all. Thus, for any rotary map tt, we have 

^<|7r(p)7r(g)|<3H. 

for all points p,q e P U Q. The lemma follows immediately. □ 

Lemma 3.8. There is a set of 0(1) rotary maps {7ri,7r2, . . . ,TTk} such that any crossing edge of 
PUQ is mapped to a relaxed crossing edge of 7rj(P U Q) by some tTj. 

Proof: Rotations about the z-axis are rotary maps, but since they do not actually change the 
radius of any sphere, we would like to ignore them. We say that two rotary maps ki and K2 are 
rotationally equivalent if ki = p o k,2 for some rotation p about the z-axis. The rotation class of a 
rotary map k, which we denote (k), is the set of maps that are rotationally equivalent to k. Since 
any rotation class {k) is uniquely identified by the point (0,0,0) in the xy-plane, the space of 
rotation classes is isomorphic to ~ S^. 

Let Bi,B2,. . . , Bjn be the smallest empty balls containing the crossing edges of P U Q. For 
each ball Bi, let Kj denote any rotary map such that Hi{Bi) is centered on the z-axis and is not 
everted, so Ki{Bi) is an empty Delaunay ball of some crossing edge of Ki{PuQ). We easily observe 
that Hi{Bi) has radius less than 3r, so the corresponding crossing edge is relaxed. Thus, for each 
crossing edge, we have a point (kj) on the sphere of rotation classes, corresponding to a rotation 
class of maps that relax that edge. 

Our key observation is that we have a lot of 'wiggle room' in choosing our relaxing maps /tj. 
Consider the ball B of radius 3r centered at the origin; this is the smallest ball containing both 
O-P and OQ- Let W be the set of rotation classes {w) such that the radius of w{B) is at most 4r 
and w{B) is not everted. TV is a circular cap of some constant angular radius 9 on the sphere of 
rotation classes, centered at (1), the rotation class of the identity map. 

For each i, the ball Hi{Bi) lies entirely inside B, so any rotation class in W transforms Hi{Bi) 
into another ball of radius at most 4r. Thus, any rotation class in the set Wj = {{w o m) \ (w) G W} 
relaxes the ith crossing edge. Wi is a circular cap of angular radius 9 on the sphere of rotation 
classes, centered at the point (kj). 

Since each of these m caps has constant angular radius, we can stab them all with a constant 
number of points. Specifically, let 11 = {(tti), {-7^2), ■ ■ , {T^k)} C be a set of A; = 0{l/9'^) points 
on the sphere of rotation classes, such that any point in is within angular distance 9 of some 
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Figure 9. For each crossing edge, there is a constant-radius cap on the sphere of rotation classes. A constant number of 
rotation classes stab all these caps. 

point in 11. (In surface reconstruction terms, IT is a 0-sample of the sphere.) Each disk Wi contains 
at least one point in IT, which imphes that each crossing edge is relaxed by some rotation class 
(tTj) G n. Finally, to satisfy the theorem, we choose an arbitrary rotary map ttj from each rotation 

class (vTj) en. □ 

It follows immediately that PU Q has 0{r^) crossing edges. 



3.4 Charging Delaunay Edges to Volume 

In the last step of our proof, we count the Delaunay edges in an arbitrary point set S by decomposing 
it into a collection of subset pairs and counting the crossing edges for each pair. 

Let S be an arbitrary set of points with diameter A, where the closest pair of points is at unit 
distance. S is contained in a cube OS of width A. We call an edge of the Delaunay triangulation 
of S short if its length is less than 5 and long otherwise. A simple packing argument implies that 
S has at most O(A^) short Delaunay edges. 

To count the long Delaunay edges, we construct a well-separated pair decomposition of S [22], 
based on a simple octtree decomposition of the bounding cube OS. (See Agarwal et al. [2] for 
a similar decomposition into subset pairs.) Our octtree has [log2 A] levels. At each level 
there are 8^ cells, each a cube of width = A/2^. Our well-separated pair decomposition H = 
{(Pi,Qi), (P2,Q2)) ■ ■ • ) {Pm-,Qm)} Contains the points in every pair of cells that are at the same 
level £ and are separated by a distance between 2>wii and Qwi. 

Every subset pair (Pi,Qi) G H is well-separated: if the pair is at level i in our decomposition, 
then for some = Q{'wi), the sets Pi and Qi lie in a pair of balls of radius separated by 
distance 2rj. Thus, by our earlier arguments, the Delaunay triangulation of Pi U Qi has at most 
0(11!^) crossing edges. 

For any points p,q £ S such that \pq\ > 5, there is a subset pair {Pi, Qi) G S such that p G Pi 
and q £ Qi. In particular, every long Delaunay edge of S" is a crossing edge between some subset pair 
in H. A straightforward counting argument immediately implies that the total number of crossing 
edges, summed over all subset pairs in S, is O(A^logA) [52]. However, not every crossing edge 
appears in the Delaunay triangulation of S. We remove the final logarithmic factor by charging 
crossing edges to volume as follows. 
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Say that a subset pah' {Pi,Qi) G H is relevant if some pair of points pi G Pi and qi G Qi are 
Delaunay neighbors in S. For each relevant pair (Pj, Qj), we define a large, close, and empty witness 
ball Bi as follows. Choose an arbitrary crossing edge piqi of PiU Qi. If Pi and are at level i 
in our decomposition, then the distance between pi and Qi is at most (6 + \/^)wi. Let /3j be the 
smallest ball with pi and Qi on its boundary and no point of 5 in its interior; the radius of Pi is at 
least 2>wi/2. Let (3[ be a ball concentric with f3i with radius smaller by wi/2. Finally, let Bi be the 
ball of radius wi inside f3[ whose center is closest to the midpoint of segment piqi. See Figure 10. 










P 
o 



o 


m 


---^ q 

o 

o 



o 




o 

o 



Figure 10. Defining the witness ball Bi for a relevant subset pair. 

Bi is clearly empty. The distance from any point in Bi to any point in S is at least Wi/2, since 
Bi C (3' . On the other hand, the triangle inequality implies that every point in Bi has distance less 
than (7 + \/3)'u;£/2 < 4.3661z/;£ either to p or to g, and thus to some cell at level (. in the octtree. It 
follows that if two witness balls overlap, their levels differ by at most [lg(7 + \/3)l = 4. Moreover, 
since any ball of radius (7 + ^/?i)wi/2 intersects only a constant number of cells at level £, at most 
a constant number of level-£ witness balls overlap at any point. 

For any relevant subset pair at level I, we can charge its 0{wf) crossing edges to its witness 
ball, which has volume Thus, the total number of relevant crossing edges is at most the 

sum of the volumes of all the witness balls. Since the witness balls have only constant overlap, the 
sum of their volumes is only a constant factor larger than the volume of their union. Finally, every 
witness ball fits inside a cube of width 8A concentric with OS. It follows that S has at most O(A^) 
long Delaunay edges. 



This completes the proof of Theorem 3.1. 



4 Extensions and Implications 

4.1 Generalizing Spread 

Unfortunately, the spread of a set of points is an extremely fragile measure. Adding a single point 
to a set can arbitrarily increase its spread, either by being too close to another point, or by being 
far away from all the other points. However, intuitively, adding a few points to a set does not 
drastically increase the complexity of its Delaunay triangulation. (In fact, adding points can make 
the Delaunay triangulation considerably simpler [26, 19].) Clearly, our results can tolerate a small 
number of outliers in the point set — up to 0(A^/n), to be precise — but this is not very satisfying. 

We can obtain a stronger result by generalizing the notion of spread. For any integer k, define 
the order-k spread of a set S to be the ratio of the diameter of S to the radius of the smallest ball 
that contains k points of S. 
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Theorem 4.1. The Delaunay triangulation of any set of points in with order-k spread Aj. lias 
complexity 0{PAl). 

Proof: The proof of Theorem 3.1 needs little modification to prove this result; in fact, the only 
required changes are in the proofs of Lemmas 3.5 and 3.6. 

Let P and Q be two sets of points contained in balls of radius r separated by distance 2r, where 
any unit ball contains at most k points. As before, we separate P and Q with a grid of 0{r'^) 
circular pixels of constant radius e. Lemma 3.4 applies verbatim. 

A simple modification of the proof of Lemma 3.5 implies that if e < 1/16, then each point of P 
is the cndpoint of at most k relaxed edges through each pixel vr. Specifically, if p is the endpoint of 
k + 1 edges pqo, pqi, . . . ,pqk that all intersect vr, then some pair of points qi and qj would be more 
than distance 1 apart, which implies that either pq^ or pqj is not relaxed. Similarly, since at most 
k unit-diameter balls overlap at any point, the set of relaxed endpoints E C P contains at most 
0{kr) points; the rest of the proof of Lemma 3.6 in unchanged. 

It follows that at most 0{k^r) relaxed edges intersect any pixel, so there are 0{k'^r^) relaxed 
crossing edges overall. Lemmas 3.7 and 3.8 now imply that there are 0{k'^r^) crossing edges, and 
the well-separated pair decomposition argument in Section 3.4 completes the proof. □ 

This generalization immediately implies the following high-probability bound for random points. 
Surprisingly, this result appears to be new. 

Corollary 4.2. Let S be a set of n points generated independently and uniformly in a cube in M^. 
The Delaunay triangulation of S has complexity 0(n log n) with high probability. 

Proof: Let C be a cube of width (n/ In n)^/^. If we generate S uniformly at random inside C, the 
expected number of points in any unit cube inside C is exactly In n, and Chernoff 's inequality [83] 
implies that every unit cube inside C contains O(logn) points with high probability. Thus, with 
high probability, = 0((n/ In n)^/^) for some k = O(logn). The result now follows immediately 
from Theorem 4.1. □ 

4.2 Unions of Several Dense Sets 

We can also generalize our upper bound to sets that do not have small spread, provided they can 
be decomposed into few subsets, where each subset has low spread. If all the subsets have the same 
'scale', the upper bound is almost immediate. 

Theorem 4.3. Let Pi, P2, . . . , P^ be sets of points in M^, each with closest pair distance at least 1 
and diameter at most A. The Delaunay triangulation of PiU P2U ■ ■ ■ U Pk has complexity 0(/c^A^). 

Proof: It suffices to consider the case k = 2; for larger values of k, we separately count Delaunay 
edges for all (2) pairwise unions. 

Let P and Q be two sets, each with closest pair distance at least 1 and diameter at most A. 
We say that an edge in the Delaunay triangulation off P U Q is bichromatic if it joins a point in 
P to a point in Q, and monochromatic otherwise. Theorem 3.1 immediately implies that there are 
O(A^) monochromatic edges. 

If P and Q are well-separated, every bichromatic edge is a crossing edge, so by our earlier 
analysis, there are O(A^) bichromatic edges. Otherwise, we can define a well-separated pair de- 
composition by building an octtree over the bounding box of PUQ, which has width at most 4A, 
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so that every bichromatic edge is a crossing edge for some well-separated subset pair. By charg- 
ing bichromatic edges to empty witness balls exactly as before, we conclude that the number of 
bichromatic edges is still O(A^). □ 

This also provides an alternate proof of Theorem 4.1, since any point set whose order-A; spread 
is Afc can be partitioned into 0{k) subsets satisfying the conditions of the theorem. 

With more effort, we can establish a similar upper bound for unions of arbitrary low-spread sets 
with arbitrarily different scales. As usual, we start by considering the case of two sets contained 
in disjoint balls. Let P be a set of points with closest pair distance 1 inside a ball QP of radius r, 
and let Q be a set of points with closest pair distance (5 S> 1 inside a ball QQ of radius R. We say 
that P and Q are well-separated if the distance between QP and QQ is at least r + R. 

First suppose P and Q are well-separated. To analyze the number of crossing edges, we follow 
precisely the same outline as our earlier proof. After an appropriate rigid motion, QP is centered 
at (0, 0, 2r) and that QQ is centered at (0,0, —2ii). We place a grid of O(r^) circular pixels of 
radius e = 0(1) on the plane z = 0, so that every crossing edge passes through a pixel. The proof 
of Lemma 3.4 immediately implies that the crossing edges passing through any pixel lie in a slab 
of width 0{R/r) between two parallel planes. 

Say that an crossing edge is relaxed if its endpoints lie on an empty sphere of radius 0{R + r). 
The proof of Lemma 3.5 immediately implies that each point in Q is an endpoint of at most one 
relaxed edge passing through any pixel. (However, a point in P might be an endpoint of more than 
one relaxed edge if R/6 < r.) Lemma 3.6 generalizes as follows. 

Lemma 4.4. The relaxed edges passing through any pixel tt are incident to at most 0{R/6 + 
R^ /rd"^) points in Q. 

Proof: Let a be the slab of width 0{R/r) containing the crossing edges through tt, and let a' be 
a parallel slab of width 0{R/r + 5) with the same central plane. We define the 'Swiss cheese holes' 
H exactly as in the proof of Lemma 3.6. For each relaxed endpoint q & Q, let Ug be a ball of radius 
6/2 centered at q; these balls are disjoint. After an appropriate scaling. Claim 3.6.1 implies that 
the surface area of H is 0{{R + r){R/r + S)) = 0{R^/r + RS), and Claim 3.6.2 implies that each 
ball Uq contains 0,(6'^) of this surface area. □ 

Since there are O(r^) pixels, there are 0{r'^R/5 + rR^/6^) relaxed crossing edges between P 
and Q. Lemmas 3.7 and 3.8 hold without modification, so the total number of crossing edges is 
0{r'^R/S + rR^/S^). 

Now suppose P and Q are not well-separated. We want to count the crossing edges — Delaunay 
edges with one endpoint in each set. Let DP be a cube of width r containing P, let C be a concentric 
cube of width 8r. Say that a crossing edge is short if both endpoints are in C and long otherwise. 
Since the spread of QDC is 0{r/S) = 0{r), our earlier well-separated pair decomposition argument 
implies that there are O(r^) short crossing edges. 

Let nQ be a cube of width R containing Q. To count the long crossing edges, we construct an 
octtree over OQ, subdividing any cell that does not lie entirely inside C, whose width is greater 
than r, and whose distance to DP is less than its width plus r. See Figure 11. Let Qi be the subset 
of Q\C inside a leaf cell £. If the width of £ is P/2', then the spread of Qi is at most R/2^5. If 
is non-empty, then P and are well-separated. 
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Figure 11. Counting long crossing edges when P and Q are not well-separated. 

By our earlier analysis, there are 0{r'^R/2^6 + rE? /A^S'^) crossing edges between P and Q^. 
Since there are a constant number of leaf cells of any particular width, there are less than 

i=l ' ^ ^ 

long crossing edges between P and Q \ C. Thus, the total number of crossing edges is 0(r^ + 
r^R/6 + rRy6^). 

The spread of P is 0(r), and the spread of Q is 0{R/5). If r < A and R/5 < A, then our upper 
bound on the number of crossing edges simplifies to O(A^). Theorem 3.1 implies that there are 
also O(A^) non-crossing Delaunay edges, so the overall complexity of the Delaunay triangulation 
of PUQ is 0(A3). 

Generalizing this analysis to more than two sets is straightforward. 

Theorem 4.5. Let Pi,P2, ... ,Pk be sets of points in R^, each with spread A. The Delaunay 
triangulation of Pi U P2U ■ ■ ■ U Pk has complexity 0{k'^A^). 

4.3 Regular Triangulations 

Regular triangulations are perhaps the most natural generalization of Delaunay triangulations. Let 
p = {p, r{p)) denote the ball centered at point p with radius r{p); we can also think of p as a point p 
with an associated weight r{p). Let S be a set of balls (or equivalently, a set of weighted points), 
and let S be the set of centers of balls in S. The power from a point x to a ball p £ S is \xp\'^ — r'^ (p) . 
The power diagram of S is the Voronoi diagram with respect to this distance function. The dual 
of the power diagram is called the regular triangulation of S. The vertices of this triangulation 
are all points in S; however, some points may not be vertices, as the corresponding region in the 
power diagram is empty. Regular triangulations can be equivalently defined as the orthogonal (or 
stereographic) projection of the lower convex hull of a set of points in one higher dimension [44]. 

The empty circumsphere criterion for Delaunay triangulations generalizes as follows. We say 
that two balls p and q are orthogonal if \pq\ = r'^{p) + r'^{q), and further than orthogonal if \pq\ > 
r^(p) +r'^{q). Any sphere that is orthogonal to a set of balls is called an orthosphere of that set. A 
subset of balls in S form a simplex in the regular triangulation of S if it has an empty orthosphere, 
that is, an orthosphere that is further than orthogonal from every other ball in S. 
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Theorem 4.6. The regular triangulation of any set of disjoint balls in whose centers have 
spread A has complexity O(A^). 

Proof: Let 5 be a set of pairwise-disjoint balls, where the minimum distance between any two 

centers is 1, and the maximum distance between any two centers is A. Note that the largest ball 
in S has radius less than A, which implies that S lies inside a ball of radius 2A. 

Say that an edge pq in the regular triangulation of S is local if \pq\ < 8min{r(p), r(g)}. We 
can charge each local edge to whichever endpoint has larger radius. By a straightforward packing 
argument, each ball p £ S is charged at most 0(r(p)^) times. The volume of each ball p is i7(r(p)^). 
Since the balls are disjoint, the number of local edges bounded by the total volume of the balls, 
which is O(A^). 

To count non-local edges, we follow the same outline as the proof of Theorem 3.1. First consider 
the well-separated case. Let P and Q be two sets of balls whose centers lie in balls of radius r 
separated by distance 2r. We claim that the regular triangulation of P U Q has 0{r'^) non-local 
crossing edges. Without loss of generality, we can assume that every ball in P U Q has radius at 
most r/4, since any larger ball has only local crossing edges. Thus, the empty orthospheres of any 
crossing edge have radius larger than 3r/4, and the bounding spheres of P and Q have radius at 
most 5r/4 and are separated by distance at least 3r/2. We modify the proof of Theorem 3.1 by 
using smallest empty orthospheres instead of empty circumspheres. This replacement increases the 
constants, but the proofs of most of the lemmas require no other modification. We will describe 
only the necessary modifications here; refer to the earlier proofs for notation and definitions. 

Edelsbrunner actually proved that regular triangulations have a consistent depth order from 
any viewpoint [42, 44]. Thus, Lemma 3.3 holds with no modification. 

The proof of Lemma 3.6 requires one qualitative change. First, we shrink each orthosphere Bp 
by only 1/8 (instead of 1/3) to obtain bp] this does not substantially affect the proof of Claim 3.6.1. 
We define a new ball Up around each endpoint p El E as follows. If r{p) < 1/2, then Up is the 
unit-diameter ball centered at p. Otherwise, Up is any ball of radius 1 /4 inside p whose center lies 
on Bp-, such a ball always exists if r > 10. A simple modification of the proof of Claim 3.6.2 implies 
that each ball Up contains 0(1) surface area of the Swiss cheese slice S. It follows that PUQ has 
0{r'^) relaxed non-local crossing edges. 

The abstract complex consisting of the regular and anti-regular^ simplices of P U Q is invariant 
under almost any conformal transformation. We can easily adapt the proof of Lemma 3.7 to show 
that for any rotary transformation tt, the spread of the centers of 7r(P) is at most a constant factor 
larger than the spread of the centers of P. (Note that the center of the transformed sphere Tr{p) 
is not necessarily the image tt{p) of the original center.) Thus, every rotary image Tr{P U Q) has 
O(r^) relaxed non-local crossing edges. Lemma 3.8 requires no modification, so P U Q has O(r^) 
non-local crossing edges. 

Finally, we slightly modify our well-separated pair decomposition argument, by defining a subset 
pair {Pi,Qi) to be relevant if and only if it contributes a non-local crossing edge to the regular 
triangulation. If this is the case, we define Pi to be the smallest empty orthosphere for any such 
non-local edge. The remainder of the argument is unchanged. We conclude that the regular 
triangulation of S has O(A^) non-local edges. □ 

We can generalize this upper bound further by allowing the balls to overlap slightly. A set of 
balls forms a k-ply system if no point in space is covered by more than k balls [81]. Applying 
precisely the same modifications as in the proof of Theorem 4.1, we obtain the following result. 

'^dual to vertices of the furthest-ball power diagram 
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Theorem 4.7. The regular triangulation of any k-ply system of balls in whose centers have 
spread A has complexity 0{k^A^). 

A special case of a A;-ply system is the hard sphere model commonly used in molecular mod- 
eling [77]. A set of balls is hard if the largest and smallest radii differ by a constant factor r 
and, after shrinking each ball by a constant factor p, no ball contains the center of any other ball. 
Halperin and Overmars [66] proved that any hard set of balls forms a k-p\y system of balls, where 
k = 0{r^p^) = 0(1). (See also Halperin and Shelton [67].) 

Corollary 4.8. The regular triangulation of any hard set of balls in whose centers have spread A 
has complexity O(A^). 

Combining the ideas in Theorems 4.1. 4.5, and 4.7, wc obtain similar bounds for any sets of 
balls that can be partitioned into a constant number of subsets, each of which is a constant-ply 
system whose centers have small constant-order spread. We omit further details. 

4.4 Surface Data 

A somewhat less obvious implication concerns dense surface data. Let S be a surface in 
M^. Recall from Section 2.2 that a point set is a uniform e-sample of S if, for some constant 
< S < 1/2, the distance between any surface point x G S to the second-closest sample point in S 
is between felfs(;r) and elfs(:r). 

Theorem 4.9. Let S be a fixed surface in M^. The Delaunay triangulation of any uniform 

e-sample ofT, has complexity 0(n'^/^). 

Proof: Let S* be a uniform e-sample of S. This set contains n = 0(/i/e^) points, where p is the 
sample measure of S [52, Lemma 3.1]. The spread of S is @{A/e), where A is the spread of S, the 
ratio between the diameter of S and its minimum local feature size. Thus, by Theorem 3.1, the 
Delaunay triangulation of S has complexity 0{A^/e^) = 0(n^/^A^//i^/^) = 0(n^/^). □ 

This bound is tight in the worst case, for example, when S is a circular cylinder with spherical 
caps [52]. Note that Theorem 4.9, which applies to any fixed surface, does not contradict our earlier 
f2(n^/ log^ n) lower bound, which requires the surface to depend on n and e. 

Attali and Boissonnat [11] recently showed that under certain sampling conditions, samples of 
polyhedral surfaces have linear-complexity Delaunay triangulations, improving earlier subquadratic 
bounds [11]. Unlike most surface-reconstruction results, their sampling conditions do not take local 
feature size into account (since otherwise samples would be infinite). They define a point set S to 
be a (e, K)-sample of S if the ball of radius e centered at any surface point contains at least 1 and 
at most K points in S; they then show that the Delaunay triangulation of any (e, K)-sample of a 
fixed polyhedral surface has complexity 0(n). 

Theorem 4.10. Let S be a fixed (not necessarily polyhedral or smooth) surface in M^. The 
Delaunay triangulation of any (e, K)-saniple ofT, has complexity 0(K^n^/^), where n is the number 
of sample points. 



Proof: Let S be an (e, K;)-sample of S. This set contains n = ri(l/£^) points, where the hidden 
constant is proportional to the surface area of E. The order-K spread of S is 0{l/£) = 0{^/n). The 
result now follows immediately from Theorem 4.1. □ 
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Again, this bound is tight (at least for constant k) in the ease of a cyhnder. Attah and Boissonnat 
conjecture that for any generic fixed surface, the Delaunay triangulation of any (e, K)-sample has 
near-hnear complexity; they define a surface to be generic if every medial ball meets the surface in 
a constant number of points. 

Finally, we consider the case of randomly distributed points on surfaces. 

Theorem 4.11. Let S be a fixed surface with surface area 1, and let S be a set of points gener- 
ated by a homogeneous Poisson process on T, with rate n. With high probabihty, the Delaunay 
triangulation of S has complexity 0{rt'/'^ log^''^ n). 

Proof: let cr be a ball of radius 0{\/ (log n)/n) centered on S. The expected number of points 

in n cj is O(logn), and Chcrnoff's inequality [83] implies that S r\ a contains between 1 and 
O(logn) points with high probability. We can cover S with less than n such balls. Thus, S is an 
(£, Ac)-sample with high probability, where e = 0{^/ (log n)/n) and k = O(logn). The result now 
follows from Theorem 4.10. □ 

Unlike most of the other results in this paper, this upper bound is almost certainly not tight. 
We conjecture that for any fixed piecewise-smooth surface, the Delaunay triangulation of a suf- 
ficiently dense uniform random sample has near-linear complexity with high probability. Recent 
experimental results of Choi and Amenta [30] support this conjecture. 

4.5 Algorithms 

Finally, our upper bounds immediately imply that several existing algorithms based on Delaunay 
triangulations are more efficient if the input point set is dense. 

Theorem 4.12. The Delaunay triangulation of any set of n points in with spread A can be 
computed in 0(A^ log ri) expected time, or in O(A^log^n) worst-case time. 

Proof: To obtain the expected time bound, we apply the standard randomized incremental algo- 
rithm of Guibas, Knuth, and Sharir [65], which inserts the points one at a time in random order; 
see also [48, 34]. The running time of this algorithm can be broken down into a point location and 
repair phases. In the point location phase, the algorithm locates the Delaunay simplex containing 
the next point. The total time for all the point location phases is 0(n log n). The repair phase 
actually inserts the point, repairs the Delaunay triangulation, and updates the point-location data 
structure. If the newly inserted point has degree k in the updated Delaunay triangulation, then its 
repair phase costs 0{k) time. 

Inserting a point into a set can only increase its spread, by either increasing the diameter, 
decreasing the closest pair distance, or both. Thus, at all stages of the algorithm, the spread of the 
points inserted so far is at most A, so every intermediate Delaunay triangulation has complexity 
O(A^). It follows that the expected degree of the ith inserted point, and thus the time for the 
ith repair phase, is 0(A^/i). Therefore, the total time for all repair phases is X^7=iO(A^/i) = 
0(A^ logn). This dominates the total point- location time. 

The worst-case bound follows immediately from the deterministic output-sensitive algorithm of 
Chan, Snoeyink, and Yap [24]. □ 
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Since the Euclidean minimum spanning tree of a set of points is a subcomplcx of the Delaunay 
triangulation, we can also compute it in 0{A'^ log n) expected time, by first computing the Delaunay 
triangulation and then running any efficient minimum spanning tree algorithm on its O(A^) edges. 
This immediately improves the 0(n^/^"^^)-time algorithm of Agarwal et al. [2] whenever A = 
0(n4/9). We can similarly improve the running times for computing other Delaunay substructures, 
such as Gabriel complexes, a-shapes [46], wrap and flow complexes [45, 57, 56, 91], and eocene 
triangles [9]. 

Theorem 4.13. Any set of n points in with spread A can he stored in a data structure of size 
O(A^logn), so that nearest neighbor queries can be answered in O(log^n) time. 

Proof: We construct a bottom-vertex triangulation of the Voronoi diagram of the points, using a 
standard randomized incremental algorithm, where the history graph is used as a point-location 
data structure. A similar algorithm for building (radial triangulations of) two-dimensional Voronoi 
diagrams is described by Mulmuley [84, Chapter 3.3]. The running time analysis is almost identical 
to the proof of Theorem 4.12. 

To speed up the search time, we add auxiliary point-location data structures to the history 
graph. Each insertion destroys several tetrahedra and creates new ones. We cluster the new 
tetrahedra according to which Voronoi cells contain them. For each cluster, we construct a point- 
location structure to determine, given a point q inside that cluster, which tetrahedron contains q. 
Because each cluster of tetrahedra shares a common vertex, we can use planar point-location 
structures with linear size and logarithmic query time [1]. With high probability, the simplex 
containing any query point q changes O(logn) times during the incremental construction of the 
Voronoi diagram. When the simplex containing q changes, we can determine which cluster to query 
by a simple distance comparison. Thus, the total time to locate q in the final Voronoi diagram is 
O(log^n) with high probability. The total space used by the auxiliary structures is bounded by 
the size of the history graph, which is 0(A^ logn) on average. In particular, some ordering of the 
points gives us a data structure of size 0(A^ logn) with worst-case query time 0(log^ n). □ 

5 Denser Regular Triangulations 

Despite our success in the previous two sections, our O(A^) upper bound does not generalize to 
arbitrary triangulations, or even arbitrary regular triangulations. Recall that regular triangulations 
in can be defined as the orthogonal projection of the lower convex hull of a set of points in R^. 

Theorem 5.1. For any n and A < n, there is a set of n points with spread A with a regular 
triangulation of complexity ri(nA). 

Proof: Any afHne transformation of R^ lifts to an essentially unique afHne transformation of R'^"'"-'^ 
that preserves vertical lines and vertical distances. Since affine transformations preserve convexity, 
it follows that any affine transformation of a regular triangulation is another regular triangulation 
(but possibly with very different weights). Thus, to prove the theorem, it suffices to construct a set 
S" of n points whose Delaunay triangulation has complexity ft{nA), such that some affine image of 
S has spread 0(A). 

Without loss of generality, assume that ^Jn|A is an integer. For each positive integer z,j < 
a/ti/A, let s(i,j) be the line segment with endpoints (2z,2j,0) ± ((-1)'+-'', (-1)*+^', 1). Let S be 
the set of n points containing A evenly spaced points on each segment s{i,j). Straightforward 



24 



Jeff Erickson 



calculations imply that the Dclaunay triangulation of S contains at least A^/4 edges between any 
segment s{i,j) and any adjacent segment s(i ib or j'lbl). Thus, the overall complexity of the 
Delaunay triangulation of S is r2(nA). Applying the linear transformation f{x,y,z) = {x,y,Az) 
results in a point set f{S) with spread 0{A). □ 

This result does not contradict Theorems 4.6 or 4.7, since the weighted points in our construction 
are equivalent to balls that overlap heavily. In fact, the largest ball in our construction actually 
contains the centers of a constant fraction of the other balls. 



6 Open Problems 

Our results suggest several open problems, the most obvious of which is to simplify our rather 
complicated proof of Theorem 3.1. The hidden constant in our upper bound is in the millions; 
the corresponding constant in the lower bound (which seems much closer to the true worst-case 
complexity) is about 10. 

We conjecture that Theorem 5.1 is tight for arbitrary triangulations. In fact, we believe that 
any complex of points, edges, and triangles, embedded in so that no triangle crosses an edge, 
has O(nA) triangles. Even the following special case is still open: What is the minimum spread 
of a set of n points in M'^ in which every pair is joined by a Delaunay edge? We optimistically 
conjecture that the answer is n/ir — o{n); this is the spread of n evenly-spaced points on a single 
turn of a helix with infinitesimal pitch. 

What is the worst-case complexity of the convex hull of a set of n points in with spread A? 
Our earlier results [52] already imply a lower bound of i7(min{A'^, nA, n^}). This bound is not 
improved by Theorem 5.1, since our construction requires points with extremely large weights. 
The only known upper bound is O(ra^). 

Another interesting open problem is to generalize the results in this paper to higher dimensions. 
Our techniques almost certainly imply an upper bound of 0{A'^) on the number of Dclaunay edges, 
improving our earlier upper bound of 0{A'^~^^). Unfortunately, this gives a very weak bound on 
the overall complexity, which we conjecture to be 0{A'^). What is needed is a technique to directly 
count [d/2j -dimensional Delaunay simplices: triangles in M^, tetrahedra in M^, and so on. 

Standard range searching techniques can be used to answer nearest neighbor queries in 
in O(logn) time using 0(ri^/ p oly log n) space, or in 0(\/n polylogn) time using 0{n) space [3, 
25, 32, 51, 75]. Using these data structures, we can compute the Euclidean spanning tree of a 
three-dimensional point set in 0(n^/^+^) time [2]. All these results ultimately rely on the simple 
observation that the Delaunay triangulation of a random sample of a point set is significantly 
less complex (in expectation) than the Delaunay triangulation of the whole set. Unfortunately, 
if we try to reanalyze these algorithms in terms of the spread, this argument falls apart — in the 
worst case, a random sample of a point set with spread A has expected spread close to A, so the 
Delaunay triangulation of the subset is not significantly simpler after all! Can random sampling be 
integrated with our distance-sensitive bounds? Is there a data structure of size 0(n) that supports 
faster nearest neighbor queries when the spread is, say, 0{^/n)? 
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